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1 Introduction 

Neutron stars are compact stellar objects with masses between around 1 and 2 solar 
masses and radii of around 10 to 15 km [I1121E]. They have magnetic helds up to 
around 10^® G [H O [HI [7j . Standard models for neutron stars traditionally assume that 
these objects are perfect spheres whose properties are described, in the framework 
of general relativity theory, by the well-known Tolman-Oppenheimer-Volkoff (TOV) 
equation m- The TOV equation is a simple first-order differential equation which 
can be solved with little numerical effort (see, for instance. Refs. [TOl fTT] j. 

The assumption of perfect spherical symmetry may not be correct. It is known 
that magnetic fields are present in all neutron stars. In particular, if the magnetic 
field is strong (up to around 10^® Gauss in the core) such as for magnetars [U 
El El Ej, and/or the pressure of the matter in the cores of neutron stars is non¬ 
isotropic, as predicted by some models of color superconducting quark matter [12], 
then deformation of neutron stars can occur [I31E11E51[I61E71E8|. We also mention 
the recent work conducted by [19] which shows that high magnetic fields in proto¬ 
quark stars modify quark star masses. The authors of this study conclude that using 
the TOV equation would be insufficient for numerical calculations of the properties 
of proto-quark stars. 
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The main goal of our study is to derive a TOV-like stellar-structure equation for 
deformed neutron stars whose mathematical form is similar to the standard TOV 
equation for spherical neutron stars. This equation will enable the user to explore the 
properties of deformed neutron stars from an equation that can be solved with rather 
little numerical effort, complementing more sophisticated numerical studies such as 
the one presented very recently in [ 20 ] , 

In contrast to the TOV stars that are composed of spherically symmetric mass 
shells, the stellar models considered in our paper are made of deformed mass shells 
which are either of oblate or prolate shape. Strategically, such a treatment is similar to 
the formalism developed by Hartle and Thorne izu, which is based on a quadrupole 
approximation of the metric of a rotating compact star. The oblate and prolate 
shapes are obtained by parametrizing the polar (z) direction of the metric in terms 
of the equatorial (r) direction along with a the deformation parameter 7 , described 
as z = yr, where we have assumed the symmetry to be axial symmetric. This 
parameter is normalized to 7 = 1 for a perfect sphere. An object that is deformed 
in the equatorial direction (oblate spheroid) is obtained for 7 < 1 , while an object 
deformed in the polar direction (prolate spheroid) corresponds to 7 > 1. Using this 
parametrization will allow us to keep the energy momentum tensor in spherical form, 
while maintaining deformation structure. 

The parametrized metric allows us to derive the stellar structure equation of de¬ 
formed neutron stars in analytic form. As already mentioned above, this equation 
constitutes a generalization of the well-known Tolman-Oppenheimer-Volkoff equation 
Pi, which describes the properties of perfect spheres in general relativity theory. 

The paper is organized as follows. In Sect.P we discuss the derivation of the stellar 
structure equation of deformed neutron stars (mathematical details are presented in 
the Appendix). The nuclear equation of state used to solve this equation is introduced 
in Sect. |3l Our equation of state is based on a relativistic nuclear lagrangian which 
describes conhned hadronic matter and a nonlocal Nambu-Jona-Lasinio lagrangian 
used to model quark deconhnement. The results are presented in Sect. HI They are 
of generic nature and do not depend on the particular choice for the nuclear equation 
of state. Conclusions are drawn in Sect. HI 

2 Stellar Structure Equations 

The properties of perfectly spherical stars are determined by the TOV equation, which 
is based on the Schwarzschild metric given by 

sm‘^{e)d(p ‘^, ( 1 ) 

where t is the time coordinate, r, 6, (p are the spatial coordinates, and <h(r) and A(r) 
denote metric functions which are determined from Einstein’s held equation of general 
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relativity theory. The deformed stellar models studied in this paper are based on a 
metric that is similar to the one of Eq. ([1]). However, instead of spherical mass 
shells the deformed stellar models are constructed from mass shells that are either of 
prolate or oblate shape. The mathematical form of the metric of such objects reads 
(G = c = l) 

ds^ = g^udx^dx'' = — dt^ + ^2 ^q2 ^2 ^ ^^ 2 '^ 

where 7 denotes a constant that determines the degree of deformation. To derive 
the hydrostatic equilibrium equation associated with Eq. (|2]), we start with Einstein’s 
held equation in the mixed tensor representation 

= -SttT^ . (3) 

Here G^i, denotes the Einstein tensor, which is given in terms of the Ricci tensor R^p, 
the Ricci scalar R and the metric tensor g^p = 6^p. The energy-momentum tensor 


T^ = (e + P) u^Up-g^pP 


(4) 


is given in terms of the stellar equation of state (pressure, P, as a function of energy 
density, e) and the matter’s four-velocity = dx^/dr and Up = dxp/dr, with the 
proper time r given by dr^ = ds^. Using Eqs. (E]) through (|1]) along with the equations 
provided in the Appendix, one arrives at the stellar structure equation of a deformed 
neutron star. 


^ _ (6 + P) [\r + dvrr^P - \r [l - 

dr (1 — 


(5) 


In the limiting case when 7 = 1 , Eq. ([5]) becomes the well-known Tolman-Oppenheimer- 
Volkoff equation [51 E] 

dP (e -I- P) (m -I- dvrPr^) 

^ = r 2 (1 _ 


which describes the structure of perfectly spherically symmetric objects. The gravi¬ 
tational mass of a deformed neutron star is given by 


dm 

dr 




(7) 


so that the total gravitational mass, M, of a deformed neutron star with an equatorial 
radius R follows as [22] 

M = 'ym{R). ( 8 ) 

In the spherical limit, the total mass is given by M = m{R) = dvry drr^e. The 
stellar radius R is dehned by the condition that pressure at the surface of a neutron 
stars vanishes, that is, P(r = P) = 0. 
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It is important to investigate the space outside the star as well. For that we 
need to examine the component of Eq. (|2]). Using the equations given in the 

Appendix, one hnds 

dr ' ’ 

One can easily see from Eq. ([9]) that asymptotically d^/dr —)• 0, as required. 

Now that we are equipped with the stellar structure equations ([5]) and ([ 8 ]), which 
are dependent on the deformation parameter 7 , we solve them for a given equation of 
state. The model chosen here assumes that neutron stars are made of quark-hybrid 
matter. It is based on a relativistic nuclear lagrangian to describe conhned hadronic 
matter and a nonlocal Nambu-Jona-Lasinio lagrangian to model quark matter [2S1 
123]. Phase equilibrium in the quark-hadron mixed phase is governed by the Gibbs 
condition. Section |3] briefly describes the key features of this equation of state. We 
stress, however, that the results presented in Sect. Hlare generic and do not depend 
on the particular choice for the nuclear equation of state. 


3 Equation of State 


3.1 Hadronic Matter 


At densities higher than that of the inner neutron star crust, and lower than re¬ 
quired for quark deconhnement, we model neutron star matter composed of baryons 
{B = {? 7 ,,p. A, S, S, A, f2}) and leptons (A = {e“,/i“}) using the relativistic mean- 
held approximation. The Lagrangian is given by uni im 125] 


C 


X] - QpT 

B 

~KmN{,gaCrf - 


■ P^) - {ruN - ga(T)\tpB + ]^{d^ad^a - 


-\ppu + - mx)^x . 

A 



( 10 ) 


The interactions between baryons are described by the exchange of scalar, vector, and 
isovector mesons {a,u,p) [26]. In the present work we employ the NL3 parametriza- 
tion as given in Tabled] [2Z]. For further details, see Refs. [101 HI] and references 
therein. 


3.2 Quark Matter 

To determine the equation of state of the deconhned quark phase we use a nonlo¬ 
cal extension of the three-havor Nambu-Jona-Lasinio model (see [23] and references 
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Saturation Properties 

NL3 Parametrization 

po (fm-^) 

0.148 

E/N (MeV) 

-16.3 

K (MeV) 

272 

m* /mjsi 

0.60 

ttsy (MeV) 

37.4 


Table 1: Parametrization of hadronic matter, where the saturation properties are 
baryonic density po) energy per baryon E/N, nuclear incompressibility K, effective 
nucleon mass and asymmetry energy a^y. 


therein). This model hosts numerous improvements over other models of deconhned 
quark matter, including but not limited to the treatment of vector interactions among 
quarks, reproduction of confinement for proper parametrization, lack of ultraviolet 
divergences with the introduction of the nonlocal form factor g{z), and momentum 
dependent dynamical quark masses. The Euclidean effective action is given by 

Se = J + m]^p{x)-^[j^{x)j^{x) - j^{x)j^{x)] ( 11 ) 

-^Tabc[ja{x)jb{x)jc{x) - 3^{x)j[{x)j^(x)] - ^jvjix)jv,fix)} , 

where 'ip = {uds)^, m = diag(m.u, m^, m^), and H, Gs, and Gy are coupling constants. 
For convenience we assume = m. The scalar, pseudo-scalar, and vector 

currents are respectively 

Jaix) = j d^zg{z)%p (^x + I) Aa^ (^~ 0 ’ 

Ja(x) = J d^zg{z)tlj + I) iysAaV' ’ (^3) 

j^x) = j d'^zg{z)tp (x + I) 

Applying standard bosonization to fllll) we derive the thermodynamic potential in the 
mean-held approximation at zero temperature |23]. We use the same parametrization 
for the nonlocal NJL model as given in Ref. [23]. The vector coupling constant {Gy) is 
given in terms of the scalar coupling constant (Gs) and is chosen to be Gy = 0.09 Gs- 
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3.3 Quark-Hadron Mixed Phase 

Phase equilibrium between the hadronic and quark phases of neutron star matter is 
governed by the Gibbs condition, 

=PQ(/^n,/ie,{V'}), (15) 

where Hn and fie are the neutron and electron chemical potentials, and { 0 } and { 0 } 
are the field variables and Fermi momenta associated with solutions of the equations of 
hadronic and quark matter, respectively. When this condition is initially met the first 
order phase transition from hadronic to quark matter begins. The relaxed condition of 
global charge neutrality allows the hadronic matter to become more isospin symmetric 
by transferring negative charge from the hadronic to the quark phase, lowering the 
asymmetry energy. This results in a mixed phase with coexisting regions of positively 
charged hadronic matter and negatively charged quark matter [ini [281 [29]. The 
equation of state for this phase is solved by combining the approaches for hadronic 
and quark matter under the Gibbs condition, baryon number conservation, and global 
electric charge neutrality. 

4 Results 

We first calculate the masses and radii of non-spherical neutron stars by solving 
Eqs. (ED and (CD numerically using the Runge Kutta method. The outcome is shown 
in Fig. [H From the results shown in this figure, we see that the maximum mass of 
the spherical (7 = 1 ) neutron star obtained for the equation of state of this work 
is 2.3 Mq. The equatorial radius of this star is close to 14 kilometers |23]. Oblate 
neutron stars are obtained for 7 values less than one, since the polar coordinate 
obeys 2; = yr. We hnd that a decrease of 7 by 10% results in a ~ 15% increase 
in gravitational mass and an increase in equatorial radius by a few kilometers (Fig. 
[ID. If we continue decreasing 7, the mass keeps increasing monotonically, ultimately 
extending into the mass region of solar-mass black holes. Prolate neutron stars are 
obtained for deformation parameters 7 > 1. In this case, as shown in Fig. [T] an 
increase of 7 by 10% leads to a ~ 12% decrease in gravitational mass and a decrease in 
equatorial radius. If one keeps decreasing 7 further, the maximum mass of a deformed 
neutron star drops down toward the 1.5 Mq region. The situation is graphically 
illustrated in Figs. [D and [5] which show the deformations of the maximum-mass 
neutron stars of Fig. [1] for 7 values ranging from 0.8 to 1.0 (oblate to spherical 
deformations) to 1.0 to 1.2 (spherical to prolate deformations). The pressure and 
energy-density profiles of the maximum-mass neutron stars of Fig. [1] are shown in 
Figs. [2] and El 

The bottom line of all this is that there may be multiple maximum-mass neutron 
stars for one and for the same model for the nuclear equation of state, depending 
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Figure 1: (Color online) Mass-radius relationships of deformed neutron stars (7 < 
1 : oblate neutron stars, 7 > 1 : prolate neutron stars, 7 = 1 : spherical neutron 
stars). The solid dots on each curve represent the maximum-mass star for each 
stellar sequence. 


on the type (oblate or prolate) of stellar deformation, which in the end is linked 
to the strengths of the magnetic helds of neutron stars and/or anisotropic pressure 
gradients in their cores. Moreover, as indicated by our calculations, the deformation 
does not need to be very large to appreciably change the bulk properties of neutron 
stars. This hnding may be critical to better understand the ever-widening range of 
observed neutron star masses and to discriminate neutron stars from solar-mass black 
holes. 

Next, we calculate the energy loss of photons emitted from the surface of a de¬ 
formed neutron star. We consider a photon created at the surface of the star (emitter) 
and leaving its gravitational held toward a detector located at inhnity, where space- 
time is hat. The photon’s frequency at the emitter, ue, is given as the inverse of the 
proper time between two wave crests, drE-, that is, ue = ^/drE = {—g^ydx^dx'')']^^'^^ 
where dx^ = dx^ = dx^ = 0 because the emitter stays at a hxed position while emit¬ 
ting the photon. The same expression written down for the receiver at inhnity reads 
^00 = l/dToo = {—g^udx^dx'^)oo ' . The ratio of these two frequencies is given by 

^00 [i-gooy^^dx^]E , . 

[{-gooyWW ^ ’ 

If we assume that the coordinate time dx^ between two wave crests is the same as at 
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Figure 2: (Color online) (a) Pressure profiles and (b) energy-density profiles in equa¬ 
torial direction for the maximum-mass stars shown Fig. [H 
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(a) 7 = 1.00; M = 2.30 M© (b) 7 = 0.90; M = 2.62 Mq (c) 7 = 0.80; M = 3.02 Mq 
Figure 4: Shapes of the oblate maximum-mass stars shown in Fig. [T] 



(a) 7 = 1.00; M = 2.30 Mq (b) 7 = 1.10; M = 2.03 M© (c) 7 = 1.20; M = 1.81 M© 
Figure 5: Shapes of the prolate maximum-mass stars shown in Fig. [T] 


the star’s surface and the receiver, which is the case if the gravitational held is static 
so that whatever the world-line of one photon is from the star to the receiver, the next 
photon follows a congruent path, merely displaced by dx^ at all points mm, this 
ratio simplihes to VooIi^e = [(— 5 'oo)^^^]e/[(— 5 'oo)^'^^]oo- Making use of the dehnition of 
the gravitational redshift, 2 ; = {ue/i^oo) — 1 , we obtain 


= 



-7/2 


- 1 . 


(17) 


Equation flTTl) shows that the gravitational redshift carries important information 
about the mass, radius, and the deformation of a neutron star. The 2 ; values at the 
equators of several deformed 1.5 M© neutron stars are shown in Fig. [1] 

The eccentricities e. 


e = sign(i?eq-^p)yl- , (18) 

of the neutron stars shown in Figs. [Hand [5] are summarized in Table [2j For spherical 
neutron stars (semi-minor axis) and (semi-major axis) are equal, so that e = 0 
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for such objects. Neutron stars whose 7 values differ by ±10% from the spherical case 
have eccentricities of e = 0.43617 if the deformation is oblate and e = —0.41601 if the 
deformation is prolate. Rapid rotation also deforms neutron stars away from spherical 
symmetry. For the neutron stars of this paper, we hnd eccentricities as low as 0.6 
for rotation at the mass shedding frequency (which sets an absolute limit on rapid 
rotation), but not smaller. The metric of a rotating neutron star can also be used 
to study the structure of deformed non-rotating neutron stars. This has been done 
recently in Ref. [20]. The results of this paper cannot be directly compared with our 
results, however, because of specihc assumptions about the energy-momentum tensor. 
We note, however, that the eccentricities of the oblate neutron stars of our study 
are compatible with those obtained in [20], depending on the degree of anisotropy 
generated by the magnetic held. 


5 Conclusions 

The goal of this work was to investigate the impact of deformation on the structure of 
non-rotating neutron stars in the framework of general relativity. For this purpose we 
hrst derived a stellar structure equation that describes deformed neutron stars. This 
equation constitutes a generalization of the well-known Tolman-Oppenheimer-Volkoff 
(TOV) equation, which describes the structure of non-rotating, perfectly spherically 
symmetric neutron stars. The mathematical structure of this generalized TOV equa¬ 
tion is such that the deformation of a neutron star (or any other compact object, 
such as a hypothetical quark star) is expressed in terms of a deformation parameter, 
7. By virtue of this parameter, models of deformed neutron stars can be built from 
non-spherical (prolate or oblate) mass shells rather than spherical mass shells. This 
leads to a stellar structure equation for deformed neutron stars which is of the same 
simple mathematical structure as the standard TOV equation and thus can be solved 
with little numerical effort. 

The parametrization introduced in our paper allows one to use a model for the 
equation of state in the limiting case of isotropy while maintaining deformation struc¬ 
ture. From the results shown in Fig. [H one sees that modest deformations can lead to 


7 = 0.80 

7 = 0.90 

7 = 1.00 

e 0.60000 

0.43617 

0 

7 = 1.00 

7 = 1.10 

7 = 1.20 

e 0 

-0.41601 

-0.55222 


Table 2: Eccentricities, e, of the oblate and prolate neutron stars shown in Figs. ID 
and (5] respectively. 
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appreciable changes in a nentron star’s gravitational mass and radins. In particnlar, 
we find that the mass of a nentron star increases with increasing oblateness, bnt de¬ 
creases with increasing prolateness. This opens np the possibility that, depending on 
the degree of stellar deformation, there may exist mnltiple maximnm-mass nentron 
stars for one and the same model for the nnclear eqnation of state, which is drastically 
different for spherically symmetric nentron stars whose mass-radins relationships are 
characterized by one and only one maximnm-mass star. This finding may be critical 
to properly nnderstand the ever widening range of observed nentron star masses and 
to discriminate nentron stars from solar-mass black holes. 
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6 Appendix 


Below, we outline the derivation of the stellar structure equation of deformed compact 
objects. For the metric given in Eq. ([2]), the non-vanishing Christoffel symbols are 

= F‘^^,J, = - , = cot{6 ), F^’^.J, = -13 rsm^{9 ), = - sin(6') cos(6'), (19) 

r 

where primes denote derivatives with respect to the radial coordinate, r, and 

/ 3 ^ , ( 20 ) 

The components of the Ricci tensor for the metric of Eq. ([2]) are calculated to be 
1 


R\ = 


/9<h'(r)m'(r)7r — <F'(r)m(r )7 — (<h'(r))^ r m{r) — $"(r)r^ 

( 21 ) 


r(r — 2m(r)) i 
-|-2<F"(r)rm(r) — 2<I>'(r)r -|- 4<I>'(r)m(r) 
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_ 

r‘^{r — 2m{r)) 


/9$"(r)r^ — 2$"(r)r^m(r) + ($'(r))^ — 2 ($'(r))^ r^m{r) 




— 7 $'(r)m'(r)r^ + 7 $'(r)rm(r) — 2'ym'{r)r + 2'ym{r) 


( 22 ) 


R\ = 


r^(r — 2m{r)) 

+ r — 2m{r) — (3^ + 2f3m{r) 


(3r‘^^\r) + 2/97r$'(r)m(r) + (3'ym\r)r — 


and 


p0_ pS 


(23) 

(24) 


The Ricci scalar, R, is calculated to be 
2 


R = 


/97<h'(r)m'(r)r^ — /37<h'(r)m(r)r — (3 (<l>'(r))^r^ 


r2(r — 2m(r)) 

+2/9 (<l>'(r))^ r^m(r) — /9$"(r)r^ + 2/9$"(r)r^m(r) — 2m{r) — 2/97^<h'(r) 
+4/9$'(r)m(r) + 2(3'ym'{r)r — 2(3'ym{r) + r — 2m(r) — /37 + 2(3m{r) . (25) 

Substituting Eqs. fET]) to fl2T)) along with Eq. fl25]l into Einstein’s held equation (E]), 
one arrives at the general relativistic stellar structure equation (j5]) of deformed com¬ 
pact objects. 
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